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ABSTRACT 

We present new reddening measurements towards the young supernova remnant Cassiopeia A, using 
two techniques not previously applied to this object. Our observations of the near-infrared [Fell] 
1.257/xm and 1.644/xm lines show the extinction to be highly variable across the remnant, increasing 
towards the west and the south, consistent with previous radio and X-ray observations. While the 
absolute value of Ay as determined by the [Fe II] lines is uncertain due to conflicting calculations and 
observations of their intrinsic flux ratio, parts of the remnant without previous optical measurements 
show comparatively higher reddening. We find Ay = 6.2 ± 0.6 from the broadband shape of the 
infrared synchrotron emission of a knot within 13" of the expansion center. Given this reddening, the 
apparent faintness of the nascent supernova, and iron mass constraints from X-ray observations, we 
estimate an ejected mass of 56 Ni of 0.058 - 0.16M©. Taken with 7-ray observations of the 44 Ti decay 
chain, this nickel mass is broadly consistent with the solar Ca/ 56 Fe ratio. 

Subject headings: supernova remnants — supernovae: individual (SN1680) — dust, extinction - 
nuclear reactions, nucleosynthesis, abundances 



1. INTRODUCTION 

Cassiopeia A is the remnant of one of the most re- 
cent known supernova in the Galaxy. It is the prototype 
of the "oxygen-rich" class of supernova remnants (OS- 
NRs), whose optical spectra are dominated by highly 
Dopplcr shifted lines of oxygen and other advanced burn- 
ing products, but are nearly devoid of hydrogen and he- 
lium. Because these remnants have not yet swept up 
significant interstellar material, their abundance pattern 
reflects that of their progenitor star and supernova. As 
such, the OSNRs, and Cas A in particular, allow unpar- 
alleled access to the physics of massive stars and super- 
nova, and should provide tight constraints on models of 
nucleosynthesis, stellar mixing, supernova physics, and 
ISM enrichment. 

Given its distance (3.4 j^'^ kpc, Reed et al. 1995), 
circumpolar location in the northern sky, and the most 
quoted reddening towards Cas A (Ay ~ 5, Hurford & 
Fesen 1996, hereafter HF96), it is surprising that its 
nascent supernova was largely unreported. SN1604 (Ke- 
pler's SN), SN1572 (Tycho's SN), SN1054 (the Crab Neb- 
ula) and SN1006 were all widely observed and recorded 
by contemporary astronomers in Europe, the Far East, 
or both. Cas A was not. While Flamsteed may have 
seen the event as a 6th magnitude optical transient in 
1680 (Ashworth 1980, though Kamper 1980 provides a 
dissenting opinion) , it was certainly less remarkable than 
other historical Galactic supernovae. Why? Was the Cas 
A supernova an unusual, under-luminous event? Were 
contemporary observers unlucky, and the supernova at 
peak brightness was highest in the northern sky during 
daylight? Is the extinction towards Cas A greater than 

1 Visiting Astronomer, Kitt Peak National Observatory. KPNO 
is operated by AURA, Inc. under contract to the National Science 
Foundation. 



previously assumed? 

Aside from this historical oddity, the apparent faint- 
ness of Cas A's outburst presents an astrophysical co- 
nundrum. The peak luminosity and characteristic light 
curve of supernovae are driven by energy injected into 
the expanding stellar debris by the radioactive decay of 
56 Ni and its daughter nucleus 56 Co. Thus, the peak 
brightness is proportional to the ejected mass of 56 Ni. 
Since the complete explosive silicon burning that pro- 
duces radioactive nickel takes place in the deepest layers 
of the supernova, the nickel yield is intimately linked with 
the explosion mechanism, energy, mass cut, and fall-back 
fraction. Its diagnostic power is further enhanced when 
compared to a supernova's 44 Ti yield. This trace isotope 
is produced by "a-rich freeze-out," a process of incom- 
plete silicon burning that dominates in regions of higher 
entropy (i.e. lower density). The high 44 Ti abundance, 
as determined by 7-radiation observations of its decay 
chain (Iyudin et al. 1994, Vink et al. 2001) relative to 
the implied 56 Ni mass inferred from Cas A's apparently 
faint outburst, are inconsistent with yields expected from 
symmetric explosion calculations (Timmes et al. 1996). 
This is generally assumed to be evidence for asymmetry 
in the supernova (e.g. Nagataki et al. 1998). There- 
fore, tighter measurements of 44 Ti/ 56 Ni should provide 
helpful constraints on advanced supernova explosion cal- 
culations. Furthermore, 44 Ti appears to be produced 
exclusively in core collapse and sub-Chandrasekhar mass 
SNIa (Timmes et al. 1996). A measurement of radioac- 
tive titanium and nickel in Cas A may inform Galac- 
tic chemical evolution and stellar population synthesis 
models in determining the relative importance of these 
two types of supernovae in producing the observed solar 
abundance ratio of 44 Ca/ 56 Fe, the two stable daughter 
products of these nuclei. 

Thus a proper inventory of 56 Ni or, equivalently 56 Fe, 
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in Cas A is desirable. Willingale et al. (2003), us- 
ing XMM-Newton X-ray observations, find an X-ray 
emitting iron mass (presumably dominated by 56 Fe) of 
0.058M©. Because the cooling time of the X-ray emitting 
plasma is generally longer than the age of the supernova 
remnant, this is, to a good approximation, the sum total 
of iron that has been heated to X-ray emitting tempera- 
tures upon passage through the reverse shock. However, 
it does not include denser material that retained enough 
electrons to cool efficiently through EUV, optical, and 
infrared emission. The cooling times for these plasmas 
can be weeks or months, making a full inventory of gas 
of these densities at reverse shock passage difficult. Sim- 
ilarly, material that has not yet encountered the reverse 
shock is nearly invisible. 

Here we take a different approach and attempt to 
bound the allowed total mass of ejected 56 Ni. Given 
that Cas A's SN was apparently unremarkable at Earth, 
knowledge of the extinction towards the supernova, a 
distance measurement, and a limit on the apparent mag- 
nitude place constraints on the intrinsic peak luminosity 
of the event. Previous optical measurents of the redden- 
ing suffer from a selection effect, in that the most com- 
mon method for estimating Ay is infeasible in regions 
of higher extinction. Existing CO (Troland et al. 1985) 
and X-ray (e.g. Keohane et al. 1996) observations show 
the absorption to Cas A to be spatially variable, with a 
general trend of increasing extinction from east to west, 
though poor spatial resolution washes out structure on 
scales smaller than the beam and will underestimate Ay 
if the absorbers are clumped. Indeed, interferometric ob- 
servations of H2CO absorption towards Cas A (Reynoso 
& Goss 2002) show considerable dumpiness. Moreover, 
radio and X-ray measurements of the absorbing column 
require an assumed gas to dust ratio to convert to optical 
extinction. As such, existing estimates of the reddening 
span several magnitudes, most commonly Ay ^4 — 8. 
(Fesen et al. 2006 assume Ay = 6 — 8 towards the com- 
pact remnant.) These uncertainties have limited the ac- 
curacy of previous attempts (e.g. Troland et al. 1985) to 
use the peak brightness to constrain the properties of the 
supernova. Here we present significantly improved mea- 
surements of Ay using [Fe II] infrared emission lines and 
a new technique based on measurement of the infrared 
spectral index of the synchrotron emission of Cas A. 

The rest of the paper is organized as follows. In §2, 
we detail the new near-infrared images and spectroscopy 
and the archival mid-infrared imaging we use for our 
measurements. In §3 we address the challenges of mea- 
suring the extinction towards Cas A with emission lines, 
estimate Ay from our [Fell] observations, and describe 
and present results from our measurement of the infrared 
synchrotron emission. In §4, we bound the ejected nickel 
mass and discuss its implications. §5 is a summary of 
our results. 

2. OBSERVATIONS 

2.1. Near- Infrared Imaging 

We imaged Cas A with the PISCES (McCarthy et al. 
2001) near- infrared camera attached to the Steward Ob- 
servatory 2.3-meter Bok telescope, over two different runs 
in summer 2002 and 2003. The 8/5 circular field of view 
and 0."5 pixels at the Bok f/9 focus allowed us to image 
the entire remnant in a single pointing while adequately 




Fig. 1. — Our color JHK S (red, green, blue) PISCES image 
of Cassiopeia A. The clumpy blue/green features are line emis- 
sion from shock heated ejecta and circumstellar material, while 
the smoother red features are synchrotron emission. Synchrotron 
"Knot N" from Wright et al. (1999) is clearly visible as a diffuse 
red feature near the center of the remnant. 

sampling typical Kitt Peak seeing, which ranged from 
0."8-l."5. All observations were made with broadband 
filters from the 2MASS system. 

In July 2002, we used a 3x3 dither pattern, with 60 sec- 
ond individual exposures in J, and 30 second exposures 
in H. The total exposure times were 54 minutes in J and 
40 minutes in H. In July 2003, we used a similar strategy, 
observing Cas A in J (45s individual, 45 minute total ex- 
posure), H (12s individual, 35 minutes total exposure) 
and K s (10s individual, 56 minutes total exposure). 

Standard image processing, including dark subtrac- 
tion, flat-fielding, distortion correction, and image stack- 
ing were accomplished with the IRAF 2 package. For 
the 2002 run, flat-field frames were constructed from 
median-filtered off-source dithers interspersed with on- 
source pointings. In 2003, no off-source pointings were 
taken, and flats were made from source frames. This re- 
sulted in a poor quality flat-field, and an underestimation 
of the flux in low-surface brightness, high-coverage fea- 
tures like the diffuse synchrotron emission observed in all 
bands. As such, the 2003 J and if images were not used 
for flux measurements, though they were used to estimate 
a correction factor for the 2003 K s data. Magnitude cali- 
bration was achieved by comparison with 2MASS (Skrut- 
skie et al. 2006) stars in the observed field. 

We present in Figure 1 a color JHK S image of Cas A 
from our PISCES data. The clumpy shocked material 
has strong lines in the J and H bands and appears blue- 
green, while the K s image is dominated by smoother 
synchrotron continuum. While this non-thermal emis- 
sion is detected in all bands, it is strongest at longer 
wavelengths, and displays red. 

2.2. Near-Infrared Spectroscopy 

2 IRAF is distributed by the National Optical Astronomy Ob- 
servatory, which is operated by the Association of Universities for 
Research in Astronomy (AURA) under cooperative agreement with 
the National Science Foundation. 
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Using our J band image to target knots of interest, 
we observed four slit positions on 26-27 September 2004 
with the FLAMINGOS 3 multi-object spectrograph on 
the KPNO 4-meter telescope. The JH grism and 6 pixel 
wide long-slit provided simultaneous wavelength cover- 
age from approximately 0.95-1.8/mi across a 10' x 1.9" 
slit on the sky, and delivered approximately 25A resolu- 
tion across the J and H bands. The emission knots in Cas 
A are densely packed on the sky, so it was impossible to 
use the standard 10" nod along the slit for sky subtrac- 
tion. Since the FLAMINGOS slit is approximately twice 
as long as the diameter of Cas A's bright ring, we placed 
the all supernova remnant line emission at one end of 
the slit, and used a 5' nod. While this large throw cost 
extra overhead (including reacquisition of a guide star at 
each nod), it ensured knots of interest would not over- 
lap between nods. Utilizing an "AABB" pattern, each 
individual exposure was 300 seconds, for a total A+B 
on-source exposure time of 4800-6600 seconds. 

Standard near-infrared spectroscopic reductions (dark 
subtraction, flat-fielding, sky subtraction, distortion 
correction) were accomplished with a combination of 
IRAF and custom Perl Data Language 4 (Glazebrook 
& Economou 1997) scripts. For the J band, we set 
the wavelength calibration with HeNeAr lamps. H band 
wavelength calibration was made against the night sky 
lines, using VLT/NIRMOS OH line list and sky spec- 
trum convolved to our resolution (Rousselot et al. 2000). 
The instrumental throughput and telluric transmission 
was measured with observations of the nearby G2V star 
HD2 12809, divided by the synthetic solar spectrum in- 
cluded with the spectral synthesis code SPECTRUM 
v2.73 (Gray & Corbally 1994). Small corrections to the 
shape of the instrumental response were corrected with 
observations of the A0V star HD240290, relative to a 
synthetic Vega spectrum calculated with the same code. 
Finally, the absolute flux scale was set with an obser- 
vation of G191B2B, utilizing the STIS/NICMOS fluxes 
(Bohlin 2007). 

2.3. Spitzer Mid-Infrared Imaging 

We used archival Spitzer/IRAC images of Cas A orig- 
inally presented by Ennis et al. (2006). While data ex- 
ists for all four IRAC bands, the broader point spread 
function for the longer wavelengths degrades our ability 
to disentangle the synchrotron emission peaks from both 
the emission-line emitting regions and the smooth contin- 
uum background, so we limit our analysis to channels 1 
& 2 (3.6//m & 4.5/im, respectively). Data products were 
downloaded from the archive, processed and mosaicked 
using standard IRAC processing. 

3. THE EXTINCTION TOWARDS CAS A 

3.1. Reddening Measurements from Emission Lines 

Reddening measurements of emission-line nebula re- 
quire observations of at least two lines of known intrinsic 

3 FLAMINGOS was designed and constructed by the IR, in- 
strumentation group (PI: R. Elston) at the University of Florida, 
Department of Astronomy, with support from NSF grant AST97- 
31180 and Kitt Peak National Observatory 

4 The Perl Data Language (PDL) has been developed by K. 
Glazebrook, J. Brinchmann, J. Cerney, C. DeForest, D. Hunt, T. 
Jenness, T. Luka, R. Schwebel, and C. Soeller and can be obtained 
from http://pdl.perl.org 



ratio and significantly different energies. With the adop- 
tion of a reddening curve (e.g. Cardelli et al. 1989, 
which we use here) and an assumption about general-to- 
selective extinction (Ry = Ay /E(B — V) = 3.1 is the 
standard value for the diffuse ISM), the observed fluxes 
of the selected lines gives the extinction. Ha and H/3 
(and higher order lines of the Balmer series) are most 
often used since they are generally among the bright- 
est lines in nebular spectra and, in the usually appli- 
cable case where their emission is primarily by recom- 
bination, the intrinsic ratio is well-known. However, 
these lines are problematic in Cas A for two reasons. 
First, the fast-moving knots that dominate Cas A's opti- 
cal spectra are composed of pure metals and are devoid 
of hydrogen lines. Second, while the quasi-stationary 
flocculi ( "QSFs" , shocked pre-SN circumstellar material) 
do emit in the Balmer series, the slow radiative shocks 
that illuminate these knots produce Ha by both recom- 
bination and collisional excitation, making the intrinsic 
Ha/H/3 (Balmer decrement) ratio uncertain. A value of 
3.0 is usually assumed for radiative shocks, but this in- 
creases to slower shock speeds as collisional excitation 
becomes more important. HF96 report measurements of 
the Balmer decrement for two QSFs in Cas A, but regard 
their implied values (Ay — 5.3,6.2±0.9) as upper limits. 

Alternatively, it is possible to measure the extinc- 
tion with metal line ratios. In the case where colli- 
sional de-excitation is negligible, the intrinsic flux ratio 
of lines originating from the same upper term depends 
only on their energies and transition probabilities. Un- 
fortunately, calculation of the Einstein A coefficients for 
the forbidden lines most often observed in optical/IR 
emission-line nebula is a difficult task, and results can 
vary by 20% or more for different calculations of the 
same transitions. Moreover, if collisional de-excitation 
is important the intrinsic flux ratio depends both on the 
electron density and theoretically-calculated effective col- 
lision strengths, which may have systematic uncertainties 
equal to or greater than those of the transition probabil- 
ities. 

The best extinction measurements for Cas A to date 
are from HF96, who examine the 1.03y«m blend and 
AA4069,4078A doublet of [S II], both of which originate 
in the 1 P (second excited) term of S + . (Collisional de- 
excitation of these transitions is unlikely to be important 
in the Cas A FMK's, as their critical densities are of or- 
der 10 6 cm -3 .) The issue of the accuracy of the necessary 
atomic data is apparent in HF96, who assume a ratio 25% 
greater than had been used in an earlier study by Searle 
(1971). Aside from systematics related to the atomic 
data, two practical considerations limit the utility of the 
[S II] diagnostic. As noted above, the lines of interest 
originate from the second excited term of S + , so, un- 
like the ubiquitous [S II] AA6717, 673lA doublet, at lower 
temperatures it is possible to have a significant S + ion- 
ization fraction with little population of the higher term, 
resulting in weak lines. More importantly, the 4070A 
lines are far into the blue, and are thus highly absorbed 
towards Cas A. Indeed, for Ay = 8, the highest derived 
extinction inferred from the radio study of Troland et 
al. (1985), A4070 > 11- HF96 only report reddenings in 
the northeastern bright ring of Cas A, where the radio- 
inferred extinction is the least. This is likely a selection 
effect, as it becomes increasingly difficult to measure the 
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faint 4070A lines in heavily absorbed regions. Since radio 
studies show the reddening towards Cas A to be patchy 
and variable across the remnant, in order to constrain 
more tightly the extinction at the explosion it is highly 
desirable to be able to measure reddenings closer to the 
expansion center or, at least, sample a more representa- 
tive range of sight lines. 

With the general availability of efficient near-infrared 
spectrographs, it is now possible to observe a wider se- 
lection of emission lines into regions of higher extinction. 
The strong 1.257/im and 1.644/im lines of [Fell] both 
originate in the a 4 D 7 / 2 level of Fe + , are bright in a range 
of emission-line nebulae, and in principle should be ideal 
reddening indicators for moderately absorbed sources. 
However, the intrinsic flux ratio of these important lines 
is a matter of some debate. Two modern calculations 
of the relevant transition probabilities (Nussbaumer & 
Storey 1988 and Quinet et al. 1996) are discrepant by 
30%. Because of the relatively small wavelength spac- 
ing of the two lines, this translates to a systematic error 
in Ay of several magnitudes. Hartigan et al. 2004 use 
the Quinet ratio (-F1.2575/-F1.644 = 1-04) in their mod- 
els, since it implies an Ay in HH111 that is consistent 
with optical Balmer decrement measurements (Gredel & 
Riepurth 1993). However, Nisini et al. (2005) find that 
neither the Quinet nor the Nussbaumer & Storey ratios 
produce produce Ay's consistent with optical data in 
HH1. In light of such theoretical uncertainty, an empir- 
ical determination is desirable. To our knowledge, there 
are no published 1.257/im/1.644/im flux ratios for rela- 
tively unabsorbed sources (e.g.s. The Cygnus Loop or 
the Orion Nebula), but Smith & Hartigan (2006, here- 
after SH06) present high S /N J and H band spectroscopy 
of P Cygni's nebula, and derive empirical Einstein A co- 
efficients from their observations, using the known op- 
tical reddening. This is, of course, dependent on their 
assumed extinction (Ay — 1.86, Lamers et al. 1983), 
itself reliant on the accuracy of B star NLTE model at- 
mospheres, being correct. However, given the current 
uncertainty in the actual intrinsic ratio, we accept the 
SH06 measurement (i*\. 257/-F1. 644 = 1-49) as the best 
available, and use it throughout this work. 

3.2. [Fe II] Reddening Measurements 

We extracted numerous spectra from our coadded ob- 
servations, seeking to maximize the signal to noise in the 
[Fell] lines of interest. In Figure 2, we mark the po- 
sitions of the 19 individual knots in which we measure 
the F1.257/F1.644 ratio to better than 20%. Of these, 
16 are fast-moving ejecta knots and 4 are QSFs, located 
both in the northern bright ring (where HF96 made their 
measurements) and the more highly reddened western 
and southern portions of the remnant. In Figure 3, we 
show the measured flux ratio against position angle in 
the remnant. While there is considerable scatter due to 
Poisson noise in the lines, there is significant real varia- 
tion on small spatial scales, as well as a general increase 
in reddening from east to west. Both these attributes 
are broadly consistent with with the radio measurements 
of Troland et al. (1985), who estimated Ay ^4 — 5 for 
much of the remnant, increasing to Ay ~ 8 in their west- 
ern most observation. However, they argue that much of 
the molecular gas along the line of sight must be clumped 
in dark clouds smaller than their 1.1' beam, so the ex- 
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Fig. 2.— PISCES H band image of Cas A. Our knots with 
successfully measured [Fell] ratios are marked with boxes, and 
synchrotron knot N is marked with a circle. The star symbol marks 
the expansion center as measured by Thorstcnsen et al. (2001). 
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Fig. 3. — The observed flux ratios for 19 [Fell] knots. The 
horizontal lines denote Ay assuming the empirical intrinsic flux 
ratio of Smith & Hartigan (2006). Adoption of the Nussbaumer & 
Storey (1988) ratio would decrease Ay by more than a magnitude. 
The hatched region shows the range of position angles sampled by 
HF96, and the range of reddenings from their optical [S II] lines. 

tinction would be expected to be greater in places, and 
vary on scales smaller than an arcminute. 

Aside from the striking variation in reddening, an over- 
all offset in the implied Ay as compared with the [S II] 
measurements from HF96 is clear. While there is no di- 
rect overlap in knots observed between the two samples, 
a number of our features are in the same area of the 
bright northern ring as the the HF96 observations. Yet, 
there is a marked difference of several magnitudes. The 
most likely explanation for this is the still quite uncertain 
atomic data for the forbidden lines of interest. A small 
portion of this error may come from the [S II] measure- 
ments. Our own calculation of the intrinsic [S II] ratio 
using more modern atomic rates (Keenan et al. 1993 
and Keenan et al. 1996) suggests Fi^m/F^jo — 0.58, 
a revision down by approximately 10% compared with 
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Fig. 4. — Measured infrared flux densities for Cas A synchrotron 
knot N. The solid line shows the best fit model for Ay = 6.2 and 
a = 0.80, while the dashed line is the unabsorbed model. 

HF96's adopted value. However, given the large separa- 
tion in wavelength of the two blends, this only accounts 
for a few tenths of a magnitude in Ay. Rather, the likely 
culprit is the large uncertainty in the intrinsic relative 
strength of the two [Fell] lines. Indeed, if we were to 
assume the theoretical value of Nussbaumer & Storey 
(1988) rather than the empirical ratio of SH06, the de- 
rived Ay would decrease by more than a magnitude, and 
would bring our least reddened values closer in line with 
HF96. This adjustment is problematic though, in that it 
would imply essentially zero extinction to P Cygni given 
SH06's observed flux ratios, in conflict with ultraviolet 
and optical observations. Clearly, an empirical measure- 
ment on a less absorbed object is required to fix this dis- 
tressing situation. Until its resolution, we caution that 
absolute measurement of Ay with these lines, for this 
SNR or any of the large number of astronomical objects 
that would benefit from NIR emission-line extinction di- 
agnostics, requires better atomic data. However, apart 
from an absolute measurement, the variation in the ob- 
served flux ratio is a clear indication of the high spatial 
variation of the reddening toward Cas A, and demon- 
strates that an accurate estimate of the extinction of Cas 
A's nascent supernova requires a measurement of Ay as 
close to the expansion center as possible. Our ability to 
reliably measure this ratio across the face of the remnant, 
even towards regions of comparitivcly high absorption, is 
a clear advantage over optical measurements. 

Finally, we note the existence of two knots, one FMK 
and one QSF, within 45" of the measured expansion 
center of Cas A. The reddenings for these knots are 
Ay = 7A~tf} and Ay = 6.0t\'^ (statistical errors only). 
These are the closest emission line features to the pro- 
jected explosion site to date with measured reddenings. 

3.3. Reddening Measurement from Infrared Synchrotron 

Emission 

Relativistic electrons accelerated in Cas A's shocks 
produce copious amounts of synchrotron emission, and 
place Cas A among the brightest sources in the centime- 
ter and millimeter sky. Cas A's X-ray spectrum has a 



hard, broken power-law component (Allen et al. 1997) 
in additional to its thermal plasma signature, and Chan- 
dra has revealed a number of discrete filaments with fea- 
tureless spectra (e.g. Hughes et al. 2000), which are 
generally assumed to be synchrotron emission (though, 
see also Laming 2001). It is then perhaps not surprising 
that synchrotron emission would be detected in the ~ 8 
decades of energy in between. Gerardy & Fesen (2001) 
were the first to present high-quality near-infrared im- 
ages of Cas A, and while their J band image was similar in 
morphology to clumpy, filamentary optical images domi- 
nated by emission lines, they noted the striking similarity 
of their K s image to radio continuum maps, and postu- 
lated they had detected infrared synchrotron emission. 
Jones et al. (2003, hereafter J03) confirmed this iden- 
tification with if-band imaging polarimetry of a small 
portion of the western shell. They also found the 2^m 
flux to be brighter than expected from simple extrapo- 
lation from the radio. Rho et al. (2003, hereafter R03) 
apparently confirm the "concave- up" nature of the Cas A 
synchrotron spectrum, though since their measurement is 
summed over the entire remnant, systematics like incom- 
plete star subtraction, contamination from line-emission, 
spectral index variation, and large-scale flat field erros 
are likely sources of systematic error. Indeed, Wright et 
al. (1999, hereafter W99) show significant variation of 
spectral indices for several knots in the millimeter-wave 
regime. 

Departures from a broadband power-law spectrum, ci- 
ther concave up or down, can be produced by a variety of 
plasma effects (e.g.s Eichler 1984, Reynolds 1998), and 
the detailed broadband shape of Cas A's non-thermal 
SED is still an open question. Of course, at optical 
through mid-infrared wavelengths, the situation is com- 
plicated by the substantial extinction along our sight line 
to the remnant. (For Ay = 4.5 and a standard reddening 
law, j44.5 Mm = 0.16 and is obviously greater for shorter 
wavelengths.) Here, we reverse the question: we assume 
the departure from a power-law is small over our limited 
range of energies, and determine Ay by fitting for the 
observed fall-off towards shorter wavelengths. 

Practical measurement of the infrared synchrotron 
spectral index in Cas A is complicated in that in J, H, 
and IRAC channel 2, most regions of bright synchrotron 
emission are also bright in emission lines. However, a 
number of isolated continuum knots are line free. For- 
tuitously, one of the brightest of these ( "knot N" from 
W99, marked with a circle in Figure 2), is only 13" from 
Cas A's expansion center (Thorstensen et al. 2001) mak- 
ing this the closest feature to the projected SN observed 
from the optical through mid-IR. Thus, a measurement 
of Ay at this location provides the strongest constraint 
on the actual reddening of Cas A's SN. 

In Figure 4, we plot the 1 — 5/^m flux densities for a 
circular aperture of radius 10" centered on knot N, with 
the background (determined from an annulus of inner 
radius 10" and outer radius 20") subtracted. For the 
ground-based JHK S bands, we convert 2MASS magni- 
tudes to Jy using the flux zero-points from the 2MASS 
documentation (Cutri et al. 2006). The error bars reflect 
the Poisson noise from the source and sky, and do not 
include any systematic error due to large-scale flat-field 
error or non-zero color terms in the flux calibration. We 
estimate these systematics to be sub-dominant. For the 
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IRAC bands, flux density is a product of the post-BCD 
processing. The Poisson errors are negligible. However, 
due to the extended wings in the Spitzer/IKAC point- 
spread function, there is substantial systematic error in 
the flux densities. We correct the fluxes from the images 
down by 2% (chl) and 5% (ch2), and assume a conser- 
vative 20% error in the absolute value 5 . 

In order to estimate Ay, we fit our data with a red- 
dened power law and used a Markov Chain Monte Carlo 
(MCMC) technique to sample the three-dimensional pa- 
rameter space of the fit. In Figure 5, we plot the likeli- 
hood contours in the Ay — a projection, marginalizing 
(i.e. integrating the likelihood function) over the normal- 
ization parameter. With flat priors on the parameters, 
there is a strong degeneracy between the extinction and 
spectral index, as higher Ay compensates for flatter spec- 
tra. This leads to only loose constraints on reddening. 
Indeed, the distribution of Ay marginalized over both a 
and k is quite broad and, without further constraints, we 
estimate Ay = 7.9 ± 1.9. 

The data and flat priors alone clearly allow non- 
physical values of the spectral index, and improbable 
reddenings. In order to limit the range of likelihood 
space available to the fit, it is clearly desirable to ap- 
ply more stringent priors to the parameters. While one 
could choose somewhat arbitrary allowed ranges on the 
parameters (top-hat priors), we choose instead to ap- 
ply a physically motivated prior on the spectral index, 
from the millimeter- wave radio observations. W09 report 
a = 0.80 for knot N. While they do not quote an error 
specific to this knot, they claim their general accuracies 
to be ±0.02 — 0.05. We ran a second set of Markov chains 
with a Gaussian prior with a a = 0.05. This shrinks the 
likelihood contours significantly, and tightly constrains 
the marginal probability distribution of the reddening 
(Figure 5). For this prior, we find Ay = 6.2 ± 0.6. 

Given the evidence for departures from a power-law in 
the broadband synchrotron spectrum of Cas A discussed 
previously, the actual uncertainty in the spectral index 
may indeed be greater than the statistical error quoted 
by W09. The concave-up spectra preferred by both J03 
and R03 flattens the power-law, and would translate to 
a lower index in our model. Inspection of our likelihood 
contours shows that this would drive our fit to higher Ay. 
(Indeed, with the flat priors, our data prefers a flatter a 
and more extinction.) One could construct more exotic 
priors that preferentially open parameter space towards 
lower spectral indices, but absent a strong theoretical 
guide of their functional form, and given the relatively 
paucity of data points, we rapidly approach art over sci- 
ence. Shorter wavelength data, whose extinction increase 
faster than the intrinsic spectrum drops, may have power 
to constrain more tightly either the spectral index or any 
significant departure from a power-law. Imaging of Cas 
A's synchrotron emission should be feasible in the red 
end of the optical with moderate-sized telescopes. Ide- 
ally, one would like a 0.8 — 2.5/xm dispersed spectrum 
of the emission, though this would require a significant 
investment of time on a large aperture facility. 

4. DISCUSSION 

4.1. Cas A's Supernova at Maximum Light 

5 http: / /ssc. spitzcr.caltech.edu/irac/calib/extcal/index. html 
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Fig. 5. — Likelihood contours for Cas A knot N for flat priors 
(red) and a Gaussian prior on the spectral index (green). The 
contours correspond to 68%, 90%, 95%, and 99% confidence for 
each set of priors. The horizontal line is the millimeter spectral 
index from Wright et al. (1999, solid), and their la errors (dashed). 

Accepting our new extinction measurement from syn- 
chrotron knot N as the most likely Ay toward Cas A's 
nascent supernova, we now seek to constrain a range of 
possible intrinsic peak luminosities, and therefore, 56 Ni 
yields. Following Young et al. (2006), 

Lpeak — MnL ©(^peafc) A(t peak ) (1) 

where Mjvi is the ejected mass of 56 Ni, Q(t) is the instan- 
taneous energy decay rate of 56 Ni and 56 Co, and A(t) is 
the efficiency with which that energy is deposited in the 
supernova gas. O(i) is defined as 



0(t) = 



Na 
56 



E 



Ni 



-t/rm 



E, 



Go 



TNi 



TCo — 1~Ni 



-t/raa 



,-t/ T N. 



(2) 

Here, Na is Avogadro's number, while E and r are the 
7-ray energy and mean decay lifetime of their respective 
nuclei (1.73 MeV and 7.6 x 10 5 s for nickel; 3.69 MeV and 
9.6 x 10 6 s for cobalt). We choose t pea k = 20.7 days (the 
rise time for the SNUb 2008ax, Pastorello et al. 2008), 
and A(t pea k) = 0.95 (the value for the most likely Cas 
A progenitor from Young et al., determined with 7-ray 
radiative transfer calculations). 

First, we estimate the minimum possible luminosity. 
Willingale et al. (2003) inventory 0.058 M© (no error 
bar reported) of X-ray emitting iron in Cas A from X- 
ray observations, and conclude that most of this material 
is ejecta. Assuming plausible isotopic ratios, the vast ma- 
jority of this material is 56 Fe, the stable daughter product 
of 56 Ni. If we allow Willingale's iron mass to be the sum 
total of 56 Ni ejected, the supernova's absolute magnitude 
would be My = —16.5. At a distance of 3.4+o^kpc, with 
Ay = 6.2 ±0.6, this translates to a peak apparent visual 
magnitude of my = 2.3 ± 0.7. (We have included a 25% 
error in the X-ray emitting iron mass.) 

Recently, Krause et al. (2008) observed the actual op- 
tical peak of Cas A's supernova through a light echo, and 
definitively classified it as a rare Type lib. They use the 
nearby SNIIb SN1993J, which had a peak My = -17.5 
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Fig. 6. — 56 Ni mass for two sets of assumptions for the supernova 
rise time and efficiency of 7-ray energy deposition. The horizontal 
dotted curve denotes the X-ray emitting iron mass, while the ver- 
tical dotted line represents the brightest plausible optical transient 
that might be missed (my = 2), and our best fit Ay + lo". The 
green solid line marks the allowed range of 56 Ni mass. 

(Richmond et al. 1994) as a template, assumed the 
maximum reddening (Ay = 8) inferred from the ra- 
dio, and derived a peak visual magnitude of my = 3.2. 
Our study indicates the extinction is considerably less, 
with concomitant increase in apparent brightness. There 
is, however, significant variation in the peak luminosity 
of the peculiar SNIIb class, with SN1996cb peaking at 
My ~ —16.3 (Qiu et al. 1999), so some care must be ap- 
plied in treating the SNIIb as a homogenous population. 

Clearly, Flamsteed's apparent observation of 3 Cas was 
not Cas A's SN at maximum light, though it may have 
been the supernova in decline. Krause et al. note that 
their light echo grew fainter by a factor of 18 in 140 days 
and that this ~ 3 magnitude fading matches a similar 
decline in lightcurve of 1993J. Morgan (2008) explore 
the observability of the putative SN1680 with new at- 
mospheric radiative transfer calculations and a detailed 
discussion of the timing within the year. We summarize 
his argument here, with our new extinction result. Given 
Flamsteed's 16 August 1680 6 m observation, and our as- 
sumed peak apparent magnitude of m ~ 1 — 3, SN1680's 
peak would have occurred in February- April. During this 
time, Cas A's position would have transitted during day- 
light, and would not have been visible at maximum alti- 
tude by several magnitudes, given any of our estimates of 
its apparent magnitude. However, with Cassiopeia's cir- 
cumpolar location in the northern sky, SN1680 would still 
have been quite high in both the morning and evening 
twilight. If Cas A's supernova were on the bright end of 
our estimates, while it may not have been visible during 
daylight like other historical supernovae, it would briefly 
have been one of the 20 brightest stars in the northern 
sky, and likely would have been reported. Conversely, 
if its peak magnitude were towards the fainter end of 
our range, it may have been just faint enough to escape 
notice. 

4.2. 56 Ni and u Ti 



Given our probable range of luminosity for the super- 
nova, we now seek to bound the total mass of ejected 
56 Ni. The detection of 0.058 M® of iron in the X-ray 
sets a hard floor. It is more difficult to derive a stringent 
upper limit, given the lack of a rigorous detection thresh- 
old. A m supernova would have garnered a great deal 
of contemporary attention; clearly Cas A's supernova did 
not. Would a 2 m or 3 m transient escape notice? Absent a 
clear answer, with the definition of the distance modulus 
and inverting equation (1) 

M v = m v - A v - n cx log 10 M Ni (3) 

we combine our ignorance into the my — Ay. In Figure 
6, we plot this difference against the calculated ejected 
mass of nickel. The lower limit on Mjvi, combined with 
our best fit Ay implies my < 2.4. This is uncomfortably 
bright, but not an entirely unreasonable magnitude for 
an object that may have gone unnoticed in the twilight. 
If, instead we assume at la increase over our best Ay 
and assume that any transient with m < 2 would merit 
at least a historical footnote, we place an upper limit 
on Mjvi < 0.16M Q . We regard this as a conservative 
estimate. Plausibly shorter rise times and larger energy 
deposition efficiencies can lower this limit by ~ 20%. 

Gamma-ray observations have detected the nuclear de- 
cay lines of 44 Ti daughter nuclei 44 Ca and 44 Sc (Iyudin 
et al. 1994, Vink et al. 2001) The latest 7-ray ob- 
servations imply an initially synthesized 44 Ti mass of 
1.6±oi x 10~ 4 AIq (Renaud et al. 2006), for a 44 Ti/ 56 Ni 
= 0.8 — 3.8 x 10~ 3 . This high mass and ratio is prob- 
lematic for spherically symmetric explosion calculations 
(Timmes et al. 1996), though in multi-dimensional asym- 
metric simulations some mix of differentially enhanced ci- 
nch, freezeout (Nagataki et al. 1998) and differential fall- 
back (Young et al. 2006) appear capable of reproducing a 
range of yields, some consistent with the observed value. 
Notably, the combination of the new INTEGRAL tita- 
nium measurements and our tighter bound on the ejected 
56 Ni mass shrink considerably the allowed abundance ra- 
tio space considered by Young et al (2006), and may pro- 
vide important constraints on advanced supernova explo- 
sion calculations. Finally, we note that our radioactive 
Ti/Ni bound is consistent with the solar 44 Ca/ 56 Fe ratio 
(1.5 x 10 -3 , Anders & Grevesse 1989) the stable isotopes 
of the considered decay chains. While is it dangerous to 
extrapolate from one observation, if this result is gener- 
ally valid, Cas A-like supernovae may lessen the need for 
the sub-Chandrasekhar mass SNIa required by Timmes 
et al. (2006) to match the solar abundance. 

Finally, we note an alternative explanation for the ap- 
parent faintness of Cas A's supernova. While nearly 
all long duration gamma-ray bursts show late-time su- 
pernova light curves, a small number do not (e.g.s. 
GRB060606 and GRB060615, Fynbo et al. 2006). Cal- 
culations by Fryer et al. (2007) show that sufficiently 
delayed explosions (after core collapse) produce a large 
mass of material on gravitationally bound ballistic tra- 
jectories, which falls back onto the compact remnant (a 
black hole for GRBs) on a timescale of minutes after the 
explosion. A correspondingly small amount of nickel es- 
capes into the remnant. In these cases, the lightcurve 
is powered not by radioactive decay, but by energy de- 
posited by the blast wave. These supernova are fainter 



and decline faster. In principle, similar such delayed su- 
pernovae could exist independently of GRBs, and their 
outbursts would be easier to hide. This idea is attrac- 
tive in that the spectra of a number of the oxygen-rich 
supernova remnants lack iron, or even oxygen-burning 
ashes such as sulfur, argon, and calcium (e.g. N132D 
and E0102, Blair et al. 2000), which may well have fallen 
back onto the compact remnant. However, given the sig- 
nificant mass of iron detected in the X-ray spectrum of 
Cas A, its supernova was unlikely to have been such an 
event, so interstellar extinction remains the likely culprit 
in its apparent faintness. 

5. CONCLUSION 

We have presented new extinction measurements for 
Cassiopeia A that show the reddening to be more vari- 
able and of higher magnitude than had been derived from 
previous optical measurements, though our estimates of 
Ay are consistent with previous radio observations. The 
[Fell] lines should be promising probes of reddening for 
regions of moderate extinction, though their use is still 
hampered by uncertainty in the intrinsic emissivity ratio 
of the two lines. Further atomic structure calculations, 
or more appropriately, observations of unabsorbed emis- 
sion line sources are necessary to yield the full utility of 



these features. Nevertheless, the observed ratio is suffi- 
cient to show the east- west reddening gradient previously 
inferred in the radio, as well as small scale point-to-point 
variations. 

We have also shown that estimation of Ay from the in- 
frared synchrotron index combined with constraints from 
the radio yields statistical errors comparable to or better 
than those from the emission lines. Fortuitously, a bright 
synchrotron knot free from emission-line contamination 
is located just 13" away from the derived expansion cen- 
ter of the remnant, and thus provides the strongest con- 
straint on Ay toward Cas A's nascent supernova. Our 
new reddening, taken with the lack of widespread re- 
portage of the supernova and with an iron mass inven- 
tory from the X-ray, bounds the ejected mass of 56 Ni 
to be 0.058 — 0.16M©. This nickel abundance and mea- 
surement of the 44 Ti yield from 7-radiation both provide 
constraints on possible asymmetry in the supernova ex- 
plosion, and imply a 44 Ca/ 56 Fe ratio consistent with so- 
lar, which may have consequences for Galactic chemical 
evolutionary models. 

Facilities: Bok (PISCES imager), Mayall (FLAMIN- 
GOS spectrograph), Spitzer (IRAC imager) 
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